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1 Introduction 

In last couple of years, amazing progress has been made by Cachazo, He and Yuan [CHY] in a series 
of papers [1-5], where tree-level amplitudes of a host of quantum field theories can be calculated using 
solutions of a set of algebraic equations. These are called the scattering eqnations and appear in the 
literature in a variety of contexts [6-14]. 

The mysterious relationship between the CHY approach and the standard QFT paradigm has been 
explained from different points of view. In [15], using the BCFW on-shell recursion relation [16, 17] the 
validity of the CHY construction for theory and Yang-Mills theories has been proven. A broader 
understanding is achieved using ambitwistor string theory [18-28], where using different world-sheet fields, 
different integrands in the CHY approach for different theories - which we will call CHY-integrands, 
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a function of the coordinates Zi in a Riemann surface - have been derived alongside with the natural 
appearance of scattering equations. A nice point of ambitwistor approach is that it provides the natural 
framework for loop scattering equations as studied in [21, 24], which lead to a breakthrough in [28]. A 
third understanding is given in [29], where inspired by the field theory limit of string theory, a dual model 
has been introduced, based on which a direct connection between the CHY approach and the standard 
Feynman diagram method has been established not only at the tree-level in [30, 31], but also at the one-loop 
level (at least for cf)^ theory) in [32] (see also [33])^. 

Although conceptually the CHY approach is remarkable and very useful for many theoretical studies 
of properties of scattering amplitudes, when applying to real evaluation, one faces the problem of solving 
scattering equations, which has (n — 3)! solutions in general. Furthermore, when n > 6, one encounters 
polynomials of degree exceeding five, rendering analytic solutions in radicals hopeless. Nevertheless, while 
the solutions can be very complicated, when putting them back into the CHY integrand and summing up, 
one obtains simple rational functions. These observations have led people to wonder if there is a better way 
to evaluate the CHY-integrand without explicitly solving the scattering equations. In [35], using classical 
formulas of Vieta, which relate the sums of roots of polynomials to the coefficients of these polynomials, 
analytic expression can be obtained without solving roots explicitly. More general algorithms are given by 
two works. In one approach [36], using known results for scalar (j)^ theory, one can iteratively decompose the 
4-regular graph determined by the corresponding CHY-integrand to building blocks related to cj)^ theory, 
thus finishing the evaluation. In another approach [30, 31], by careful analysis of pole structures, the authors 
wrote down a mapping rule, so that from the related CHY-integrand, one can read out contributions of 
corresponding Feynman diagrams. 

Both approaches are powerful and have avoided the need of solving the scattering equations explicitly. 
Furthermore, based on these perspectives, especially the mapping rule, one can use Feynman diagrams 
to construct the CHY-integrand. These results produce a very interesting phenomenon: two different 
CHY-integrands can produce the same result. For example, there are two very different CHY-integrands 
for scalar cj)^ theory: one is given in [5], while another one is given in [30, 31]. We are naturally led to 
wonder how to explain the equivalence of different CHY-integrands. 

In fact, as a rational function of coordinates Zi on a Riemann surface, the equivalence can occur on 
three different levels. 


1. At the first level, their equivalence is pure algebraic, i.e., through some algebraic manipulations, one 
rational expression can be transformed to another one. For example, for 4-point amplitudes of (f>^ 


theory, on the one hand we have the integrand Ii = 


^12^23234241 212223234241 


where we have defined 


Zij = Zi — Zj which gives a contribution of ^ On the other hand, we have the integrands 

hi = 


S41 


212223234241 212224243231 


which gives-and the integrand I 22 = 


212223234241 213232224241 


which 


gives —It is easy to check algebraically that R = —hi — h 2 - Equivalences at this level is 


'^Recently, inspired by the development of CHY-approach, a new method to construct all loop integrands for general 
massless quantum held theories has been proposed in [34]. 
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of course rather trivial and in order to proceed to the other two levels of equivalences, we need to 
change our viewpoint to algebraic geometry, i.e., to transform the scattering equations to a set of 
polynomials of (n — 3) variables, defining an ideals X; 

2. The difference of two CHY-integrands can be written as J(/i — 12 ) = ^ where both P, Q are polyno¬ 
mials and J is the Jacobian we will review shortly. If P belongs to the ideal X, then for each solution 
of the scattering equations P = 0, thus at the second level we say that Ii is equivalent to 12 ', 

3. However, in practice, most of the time something more complicated happens and we find that though 
P does not belong to the ideal X and J{Ii — I 2 ) = 0 when and only when we sum over all solutions. 
If this happens, we say that Ii is equivalent to I 2 at the third level. It is clear that this is the most 
involved situation, and indeed, in practice this is the most frequently encountered. 

Motivated by the above considerations and bearing in mind that indeed the most conducive perspective 
on studying the scattering equations is through the language of algebraic varieties and polynomial ideals 
[37, 38], we turn to this method of attack. The above problem thus translates to finding the sum over 
the rational function ^ evaluated at the roots of a zero-dimensional ideal I, and testing whether the sum 
is zero. Luckily, there is a theorem in commutative algebra, due to Stickelberger, which addresses the 
situation [39]. We will discuss the theorem and the associated algorithm in illustrative detail. It turns 
out that this method not only checks the equivalence at the third level, but also evaluates the integration 
without solving the scattering equations. In this sense, it is in the spirit of the methods in [36] and [30, 31]. 
Although it is sometimes less efficient compared to these two methods, it does provide a very different 
angle to approach the problem and could have very advantageous repercussions. 

The structure of the paper is as follows. We begin with a brief review of the tree-level scattering 
equations in §2, before laying down the foundations of the theory of zero-dimensional ideals in §3, especially 
that of companion matrices. We then illustrate the technique with ample computational examples in §4, 
before concluding with remarks in §5. 

2 Review of Tree-level Scattering Equations 

In this section, we offer a brief review of tree-level scattering equations and the reader is referred to [1-3, 5] 
for details. The scattering equations are given by 

X, = = a = l,2,...,n, (2.1) 

tf Za- Zb 
b^a 

where Sah = {ka + = ‘^ka ■ kb-, and ka with a = 1, 2,..., n are n massless momenta for n-external particles 

and Zi are complex variables living on with n punctures. Although there are n equations, only (n — 3) 
of them are linear independent after using the momentum conservation and massless conditions which 
translate to the following three relations 

^Xa = 0, '^£aZa = Q, ^ = 0 , 
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a 


a 
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which are, in fact, the consequence of the SL{2,C) symmetry on the CP^. Because of this, we can insert 
only (n — 3) delta-function. To make sure the result does not depend on which three equations have been 
removed, we make following combination and define the measure^ 


— Zij5 (£a) 1 

a^i,j,k 


(2.3) 


with Zij = Zi — Zj. With the above, the general tree-level amplitude is given by 




n 


nr=i dz. 

vol(5L(2,C)) 


Q{S)T{z) 




(2.4) 


where duo = '^^rdzsdzt ^Qj^^gg after we use the Mobius SL(2, C) symmetry to fix the location of three of 
the variables Zr,Zs,zt by the Faddeev-Popov method. Different QFTs give different forms of the CHY- 
integrand J^{z). Invariance under the Mobius transformation requires J-{z) to have proper transformation 
behaviors, i.e., under z' = we have 


F{z) 


/-r-r {cZj -b d)^ \ 

[l{{ad-bcr) 


F{z) . 


(2.5) 


To simplify expression (2.4) further, we integrate out the delta-functions to arrive at the key expression 


sol 


^ij^jk^ki^rs^st^tr ^ 
2+j+Zc+r+s+t I $ I ’ 


(-) 


( 2 . 6 ) 


where three arbitrary indices i,j, k correspond to three removed scattering equations while three arbitrary 
indices r,s,t correspond to the above mentioned three fixed locations. The sum is over the solution set 
of the scattering equations, which is generically a discrete set of points. Furthermore, in the above, the 
Jacobian matrix is calculated as (a for rows and b for column) 


^ab 


dzb 


^ab 



a ^ b 
a = b ’ 


(2.7) 


and is the determinant of after removing the i-th, j-th and k-th rows and r-th, s-th and t-th 

columns. 


Specific Examples: Now we list some examples in the literatures [2, 3] (more can be found in [5]). 

According to the CHY formula, the integrand unifying scalars(6 = 0), gluons(6 = 1) and gravitons (6 = 2) 
is given by 


^f>,n — i 


Tr(r 


o(l) 


T<^{n)) . 


2-b / \b 

Pf ' ^ 


aes„/z„ ■ ■ ■ Mn)a{l) 


nice explanation of this fact can be found in [15, 37] 


( 2 . 8 ) 
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where the sum is over permutations on n elements by the symmetric group Sn, up to cyclic ordering of Z^,, 

/ _(jt \ 

is a 2n X 2n antisymmetric matrix defined by ^' = ( ^ ^ ) (where t is the transpose of the matrix), 


with A, B, C being n x n matrices with components 


C B 


{ ka'kf} 

1 

f ^a’^b 


1 Za-Z^ 

, Bab = \ 

' Za-Zb 

II 

1 0 

[ 0 


(-a-kb 

■pOT* 

^ a = h 


a ^ b 


E eg-Zi 
CT^a Za-z^_ 

and Pf ' 'I' is the reduced Pfaffian (square-root of the determinant) of 'I' defined by 

Pf ' ^ = , 

’ 


(2.9) 


( 2 . 10 ) 


where 1 <i,j <n and is the matrix removing rows i,j and columns i,j. We recall that the Pfaffian 
of a 2n X 2n antisymmetric matrix can be computed as 


= i E sgn(a)n^ 


2^n\ 


0'2i-10'2i ) 


( 2 . 11 ) 


0'GS2n 


i=l 


where sgn(cr) is the signature of cr G S 2 n- Importantly, Pf is non-zero on the solutions of scattering 
equations, while Pf ' 'I' is independent of the choice of i,j. 

Specifically, we have that 


1 


For color-ordered bi-adjoint scalar theory, 

^12^23 ^nl 

For color-ordered Yang-Mills theory with ordering {1, 2,..., n}, 

1 


-Fym = 


-Pf' . 


^ 12^23 ■■■Znl 


For gravity. 


To = (Pf ' ^)" = 


{Zi - ZjY 


Det(^^:^) . 


( 2 . 12 ) 


(2.13) 


(2.14) 


Having presented the above examples, let us go back to (2.6). As is clear from the expression, the 
right hand side is a rational function in the complex variables Zi. To employ methods developed in 
algebraic geometry, we need to associate solutions to a zero-dimensional algebraic variety defined by some 
polynomials. In other words, we should rewrite £a defined in (2.1) to an equivalent polynomial system. This 
has been done in a beautiful paper [37], where it has been shown that scattering equations are equivalent 
to following set of polynomials 

0 = hm = ) 2 < m < n — 2 , (2-15) 

SSj4,|S|=m 
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where the sum is over all subsets 5 of ^ = {1, 2,n} with exactly m elements and ks = Ylbes 

and zs = algebraic geometry, notably the affine Calabi-Yau properties of (2.15), has been 

investigated in [38]. 

A very useful observation made in [15, 37] is that If all / 0, then values of Za are all distinct. The 
set (2.15) has not fixed gauge. One of the choice of gauge will be to set, as is standard with points on CP^, 
the three points zi = oo, Z 2 = 1 and Zn = 0. Under this choice, the set of polynomial is reduced to 

hl<m<n-3^ lim ^ (ks + ki)^Zs\z2=l,z„=0 , (2.16) 

- - Zi^oo 2;i 

S&A/{l,n},\S\=m 

In summary, h defines a zero-dimensional ideal in the polynomial ring in n — 3 variables. Then, using 
the standard Bezout’s theorem, the number of points in this ideal (solutions of the scattering equation) 
is deg(^m) = (n - 3)!. 

Instead of computing the amplitude with formula (2.6) by summing over all solutions of scattering 
equations, we will show in next section that, using the companion-matrix method, we can compute the 
amplitude An = Ylsoi § ^s the trace of certain matrix composed of so-called companion matrices Tz^ 

An = Tr(F'U^T.^Q'-^U^T.J , (2.17) 

without the explicit solutions of scattering equations. 

3 The Mathematical Framework 

As mentioned in the introduction, it is expedient to consider the problem within the framework of ideal 
theory. Our problem is thus the following. 

Problem: Let I = (/*) be a zero-dimensional ideal in R = C[xi,...,Xn] generated by 

fi=i, 2 ,...,k{xi, ..., Xn) G R and let r(xi,..., x„) be an arbitrary rational function in the fraction 
field of R. Because dime / = 0, / = is a discrete set of, say N, points. We wish to 

evaluate 

N 

1=1 

where each summand is an evaluation of p at one of the discrete set of zeros Zj . In particular we 
wish to test whether this sum is 0. This is the level 3 equivalence mentioned in the introduction. 

Of course, the idea is to solve this without explicitly finding the roots Zj. This can be done using the 
technique of companion matrices [40] (cf. also [41]). Suppose a Grobner basis for I has been found for 
some appropriate monomial ordering and B is an associated monomial basis for I, which can be seen as a 
vector space of dimension d. Then the multiplication map by the coordinate variable Xi 

R/I —^ R/I 

/ -^Xif (3.1) 
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is an endomorphism of quotient rings. In the basis B of monomials, this is a d x d matrix and is called a 
companion matrix. Clearly, {Tj} all mutually commute and thus can be simultaneously diagonalized. 
We have the following [39]: 

THEOREM 3.0 (Stickelberger) The complex roots Zi of I are the vectors of simultaneous eigenvalues 
of the companion matrices i.e., the corresponding zero dimensional variety consists of the points: 

v(/) = {(Ai,..., A„) E C" : G C"Vi : TiV = A^u} . 

We point out that the original statement of the theorem is in terms of annihilators in algebraic number 
theory and is perhaps a little abstruse. Fortunately, the computational algebraic-geometry community has 
rephrased this into the readily usable form of companion matrices [40, 42]. In particular, we have the 
following important consequence: 

COROLLARY 1 Our desired quantity 

N 

Y^r{zj) = TT[r{T,,...,Tr,)] 

where the evaluation of the rational function r on the matrices Ti is without ambiguity since they mutually 
commute. 

We remark that because r is rational, whenever the companion matrices appear in the denominator, they 
are to be understood as the inverse matrix. 


3.1 Warmup 

Before proceeding to examples in our context, we present two simple exercises to demonstrate our algorithm. 
Computations can be made in Macaulay2 [42] or Singular [43], or the latter’s interface with Mathematica 
[44]. Let 

I := {xy — z,yz — X, zx — y) C R = C[x, y, z] . (3-2) 

We know, of course, that there are 5 roots 


V(/) = {(0,0,0),(-1,-1,!), (-1,1,-!), (1,-1,-1), (1,1,1)} 

Now we consider two functions, where one is polynomial and another, rational: 

3x^y + xyz 


pix,y,z) = 3x-^y + xyz, Q{x,y,z) = 


2xy‘^ + 4^2 -p 1 


(3.3) 


(3.4) 


It is easy to hnd, after summing over the solutions, that 




V{1) 


V{1) 


20 
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We now show how the companion matrices work without finding the roots (3.3) explicitly. 

In the lex ordering of x -< y -< z, the Grobner basis and the monomial basis are, respectively, 

GB{I) = {z^ - z, yz^ - y, X -yz) ] B = {1, y, yz, z, z^} . 

Therefore, we have that, in the quotient ring R/I, 

x.B = {yz, z, z‘^,y, yz] , y.B = {y, z, yz, y] , z.B = {z, yz, y, z} , 


so that 


n = 



T, = 



T. = 



(3.6) 


(3.7) 


(3.8) 


Therefore, the sum over the roots of p is 


Tr {3T^Ty + T^Tyn) = 4 


(3.9) 


and we have nice agreement with (3.4). 

For the Q, the numerator is = 3T^Ty + T^TyTz = 


/o 0 0 3 1 \ 
' 0 1 3 0 0 ' 
0 3 10 0 
0 0 0 1 3 
\0 0 0 3 1 / 


while the denominator is D = 


2T:,T^ + 4T^ + I = 


/ 1 0 2 0 4 \ 
' 0 5 0 2 0 ' 
0 0 5 0 2 
0 2 0 5 0 
\ 0 0 2 0 5 / 


Thus we calculate Tr(A^ZI ^) = which is the right answer on 


comparing with (3.4). 

Before going to examples of scattering equations, let us give some remarks. First, the theorem in its 
original form is for polynomial test functions r, while functions we will meet in scattering equations are 
rational functions, i.e., the form ^ with both P,Q are polynomials. Luckily, the theorem and corollary 
can be generalized trivially since we can diagonalize companion matrices simultaneously because the next 
remark. 

Now, there is a second part of the theorem which states that the companion matrices can be simul¬ 
taneously diagonalized if and only if the ideal / is a radical ideal. That is, there are no multiple roots. 
However, as shown in [15], if all kg 7^ 0, the solutions of Zi will all be different, so we indeed have a radical 
ideal and find simultaneous eigenvalues readily. 

Third, since there are (re — 3)! solutions, the size of Ti will be in general d = (re — 3)! which will become 
very large with re. Although with this counting, the efficiency of the method may be arguable, it does 
make the following property manifest: after summing over all solutions, the final result must be rational 
functions of k,e. 






4 Illustrative Examples 


In the following, we will use several examples to demonstrate the companion matrix method. The n = 4 
case is simple. The companion matrix is 1-dimensional, equaling to the single solution of scattering 
equations. We compute the amplitudes in scalar (j)^, Yang-Mills and gravity theories to show the validity of 
the method. For n = 5, we first study the amplitude of scalar (f)^ theory, and show that the amplitude-level 
identity can be understood by the fact that the trace of matrix is a linear mapping, and use it the explain 
a 7-point identity proposed in [31]. For the amplitudes of Yang-Mills and gravity theories, we will show 
that the companion matrix method indeed produce the correct amplitudes. 

For n = 6 , the scalar (j)^ theory will be shown to detect the pole structures so that the amplitude can 
be constructed by setting appropriate kinematics. The Next-MHV gluon amplitude is also presented as 
an example to show the validation of companion matrix method in a more difficult situation. Finally, for 
n = 7 amplitudes of scalar (jy^ theory, we demonstrate that, when companion matrices are computed in the 
diagonal form, the diagonal elements of the integrand matrix (which we recall to be an (n — 3)! x (n — 3)! 
matrix for n-points) have one-to-one mapping to the integrand computed at the (n — 3)! solutions of the 
scattering equations, so they are not only equivalent at the amplitude level, but also at the level of each 
solution as indicated by Stickelberger’s theorem. 


4.1 Four-Point Amplitudes 

The n = 4 case is trivial. There is only 4 — 3 = 1 variable left, so the companion matrix is just a complex 
number. Let us remove three scattering equations £i, £ 2 , T 4 and gauge-hx three points = 00, Z 2 = 1 and 
2:4 = 0. The remaining one scattering equation is 

536 'S13 , 523 , S34 (s 23 + 'S34)^;3 “ 534 


f3 = E 


+ 


+ 




Zy, - Zb Z3- Zi Z3 - Z2 Z3- Z4 


Z3{Z3 - 1 ) 


(4.1) 


We can define the ideal I = ((^23 + 534)23 — S34) in C[ 23 ]. It is a linear function, so the Grobner basis 
and monomial basis are trivially 

GB(/) = ((s 23 + 534)^3 - S34) , B = {1} . (4.2) 

The polynomial reduction of 23 B = { 23 } over Grobner basis of ideal I gives the remainder Thus 

in the quotient ring, the companion matrix is given by 


S 34 


= {- 

S 23 + 534 

We now proceed to the three cases of concern. 


} —t — 


S 34 


523 + 534 


(4.3) 


4.1.1 Scalar Theory 

For the 4-point amplitude in scalar theory, we wish to compute (recall that the three points 21,22,24 


have been gauge fixed) 

~2 ^2 ^2 
^ ^ 12 ^ 24^41 


sol 


|*|124 

I^ll24 


~ S3,iz3 - 1)2 + 523^1 “ Qizs) 
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(4.4) 













where we have used the simplification 


I (j) 1124 _ jjj _ '^34 

I^ll24 — ®33 —- n - 


•523 


1 


2^3 (^3 - 1 )" 


(4.5) 


^3 (^3-1)^ ^|(S 12 + S 23 ) - 2^3512 + Si 2 

SO that the factor l/ 2 ; 23^34 cancels the numerator of 1 /| 4 >|J 24 - We see that the final expression is summed 
over the (discrete) solution set of the scattering equation which is rather trivial here. The summand is a 
rational function in the free variable zs which we define as P/Q] of course, P = 1 here and Q will be used 
later. 

Finally, using the simple expression for the companion matrix from (4.3), we have 

Tr(g^) =W(- r ( ^ \ 2T -—) = = , (4.6) 

^Qi-Lzs)^ ^ - 1234 , 23(512 + S 23 ) — 2 J^ 3 Si 2 + S 12 2 S23'S34 ^14 ^12 

after some identities between Mandelstam variables have been used. This is indeed the same answer as the 
standard known result as given in the introduction. 

4.1.2 Yang-Mills Theory 

For 4-point amplitude in Yang-Mills theory, we want to compute (under gauge-fixing zi = oo, Z 2 = 1, Z 4 , = 

0 ), 

^^12^24^41 Pf'4'8x8 _ P'^^izs) 

sol 


Ar = ^ 


|<T>|124 

I^ll24 


2:i2^;232:34^;41 


So, ' 


(4.7) 


To avoid the divergence when taking the limit zi —>• 00 , one of the removed rows (columns) in T should be 
1, otherwise some terms in Pf ^Tgxs would lead to infinity. Let us then choose the reduced Pfaffian as 

-2 


Pf'4-8x8 = 


-Pf Tji 


(4.8) 


Zl - Z2 

The large zi dependence of Pf'Tgxs is then l/zf, and together with the factor from the scalar part, we 
obtain a finite integrand when taking the zi 
matrix. 


00 limit. Explicitly, the new matrix T^THisaOxG 


/ 


T = 


V 

whose Pfaffian is given by 
,12 


0 

kgki 

£lfc .3 

£ 2^,3 

CO 

£ 4^3 

Z3-Z4 

21-213 

2:2-213 

X^c^3 23-2:c 

24-23 

k^ki 

0 

£ 1^4 

£ 2^4 

£. 3^4 

\ ■' ^ 4^0 

Z4-Z3 

21-24 

22-24 

(4 

CO 

1 

i 4 

4 -^c ^4 24 — 2, 

eifca 

eifcr 

0 

eie2 

eie3 

€164 

CO 

1 

14 

2:1-2:4 

21-22 

21-23 

21-24 

€2 fes 

£ 2^4 

ei€2 

0 

£ 2^3 

€264 

2:2 -^3 

2:2-2:4 

^2-^1 

22-23 

22-24 

' £.3 fee 

£. 3^4 

eies 

^ 2^3 

0 

€364 

'C^3 Z 3 -Zc 

2:3-2:4 

23-21 

23-22 

23-24 

£ 4^3 

£ 4 fcc 

eie 4 

£2 £4 

£364 

0 

2:4--23 

4--/C^4 2:4-^C 

24-2:1 

24-22 

24-23 


\ 


(4.9) 


Pf 4 'i 2 = 'ki64'254'34 — Ti 5T264'34 “ 4'i6'k244'35 -|- 4'i4T26'k35 -|- 4'i54'24'k36 
-4'i4'k254'36 + 4'i6T23'k45 “ 4'i34'26'k45 + 'I'i24'36T45 — 'ki54'234'46 
+ 4'i34'254'46 — 4'i2T354'46 -|- 4'i 44'234'56 “ 4'i34'24'k56 + 'I'i24'344'56 . 


(4.10) 
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The reduced Pfaffian Pf in this case is a rational function with denominator z^(z 3 — l). Together with 

the factor l/ 2 ; 23^34 = — 1 ), they cancel the numerator of 1/|^>|{245 leaving a 2:3 in the denominator 

of integrand. 

Therefore, it is immediate that the numerator of the integrand comes entirely from the numerator of 
the reduced Pfaffian: 


pYM _ ^2^ _ 5^^261,362,4 — 263,4 Ki, 3K2,4 + 262,4«^1,4^3,2 — 26i,4K2,4«^3,2 + 262,4«^1,3K3,4 
+ 262,4^1,4^3,4 — 26i, 4K2,4/^3,4 + 262,3^1,3^4,2 “ 26i,3K2,3K4,2 + 261,2^3,2^4,2 
+ 26i, 2K3,4K4,2 + 26i,3K2,4«^4,3) + ^3 ( — Sl26l,462,3 + ■Sl26l,362,4 + Sl 26 l,263,4 
— 263,4Ki, 4«;2,3 + 263,4«:i,3K2,4 “ 262,4Kl,3«:3,4 — 262,4Kl,4K3,4 + 26i,4«:2,3«^3,4 


+ 26i,4K2,4K3,4 — 261,2^3,4^4,2 + 262,3^1,3«:4,3 + 262,3^1,4^4,3 “ 261,3^2,3^4,3 
— 26i,3K2,4«4,3 + 261,2^3,2«^4,3) “ Sl 26 l,263,4 , 
where 6 ^,^ = 6 j 6 j, Kij = eikj. The denominator of the integrand, on the other hand, is 

^ ^ 3 ( 5^2 + S 23 ) - 2 z|si 2 + 2:3512 = 2;3g , 
where Q is the denominator of integrand for the scalar theory from (4.4). 


In summary, by computing Tr(Pg 


-1 


we arrive at 


/(YM ~ ~ ^ ^ ^ ^ 'Sl2~ ~ 523— ~ 

-^4 — £1,362,4 — 61,462,3 — 61,263,4-61,462,3-61,263,4 

•523 S 12 


H- ( — 263 4 Ki 4 «; 2,3 + 263 4K1 3 ^ 2,4 “ 262,4^1 3K3 4 — 262,4^1,4^3 4 + 26i 4 K 2 , 3 «^ 3,4 

S12 V 

+ 261,4^2,4^3,4 — 26i, 2K3,4K4,2 + 262,3«^1,3K4,3 + 262,3Kl,4«:4,3 “ 261,3^2,3^4,3 

— 26i,3K2,4£v 4,3 + 26i,2K3,2K4,3) 

H- ( — 263,4Ki, 4K2,3 + 262,4'«1,4«^3,2 “ 26i,4K2,4'«3,2 + 26i,4K2,3K3,4 + 262,3«^1,3«^4,2 

•523 ^ 

— 26i, 3K2,3£24,2 + 26i,2K3,2K4,2 + 262,3^21,3^4,3 + 262,3/24,4^24,3 — 26i,3/C2,3/24,3 

+ 261 , 2 / 23 , 2 / 24 , 3 ^ ■ 


(4.11) 

(4.12) 


(4.13) 


The pole structures are similar to the scalar (j)^ theory, while the terms without poles come from the 
gluon four-vertex. Of course, by momentum conservation and the property €iki = 0, we can further 
write the above result as a function of all independent kinematics, for example by using identities Cjk^ = 
—ejks — ejk 2 — ejki and 64/04 = 0. This result agrees with the one computed directly by Feynman diagrams. 


4.1.3 Gravity 


For the 4-point amplitude in gravity, we want to compute 


sol 


^2 2 2 
^12^24^41 


4> 


124 

124 


Det'('I/ 8 x 8 ) 


E 


^12^24^41 


|4> 


124 

124 


(Pf'Tsxs)'^ 

Z3GSOI 


P°{Z3) 

g°(23) 


(4.14) 
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under the gauge-fixing zi = cx), Z 2 = 1, 2:4 = 0. As in Yang-Mills theory, we choose the reduced Pfafhan as 

Pf'^-Sxs = , (4.15) 

Zl — Z 2 

and as above, we know that the squared reduced Pfafhan (Pf ^'I'sxs)^ is a rational function with denomi¬ 
nator z^{z 3 — 1)^. This cancels the numerator of l/|d>|} 24 , leaving a zl in the denominator of integrand, so 
that the numerator of integrand equals to the square of numerator of reduced Pfafhan: 

pG ^ ^pYM-^2 ^ 

while the denominator of integrand is 

= 4('Si2 + S 23 ) - 22;|si 2 -h z^su = z^Q = z^Q'^^ . (4.17) 


Combining all together, we have 

pG ^pYM\2 




^pYM^2 ^pYM'^2 


(4.18) 


zsQ^^ Z 3 (QY^)2 ’ 

where Q is the denominator of integrand for scalar theory from (4.4) and the expressions for P^m 
and Q™ are given in (4.11) and (4.12). In the present case of n = 4, there is only one solution for 
scattering equations, and the companion matrix is really 1-dimensional in (4.3), so although in general 
Tr(MiM 2 ) 7 ^ Tr(Mi) Tr(M 2 ), here we simply have 


Tr 


(p) =T..(Q)T.. 




Q' 


■§12523 
512 + 523 


(M 


ym\2 


By BCJ relation [45], we can rewrite this to the familiar one 

A 2 = 5i2^r(l, 2,3,4)Mr(l, 2,4,3) , 

in agreement with the known result by KLT relation [46-50]. 


(4.19) 


(4.20) 


4.2 Five-Point Amplitudes 

For n = 5 amplitudes, there are hve scattering equations, but only two of them are independent. Under 
the gauge-hxing zi = 00 , 2:2 = Ij -^5 = 0, the Dolan-Goddard’s formula [15] gives: 

/l = 512 + 5132:3 514^4 , /2 = 5452:3 + 53524 + 5252324 . (4.21) 


We can solve these two equations to get two solutions: 


and 


— 512525 — 513.S35 + S14S45 — y/A 

soh : 23 = --- , 2:4 = 

2513 S 25 


, —512525 — 513.S35 -I- S14S45 -I- a/A 

S 0 I 2 : 23 = --- , 2:4 = 

2513 S 25 


— 512525 + 513.S35 — S14545 + 
2514S25 


— 512525 + 513.S35 — 514^45 — a/A 

2514 S 25 
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where A = (S12S25 + S13S35 — 514545)^ — 4S12S13S25S35. We can see that, in general the solutions are not 
rational functions, as is to be expected from high degree polynomials, though of course the final result 
of the sum over these points will be. One can see that the cancelations and simplifications will be very 
involved. 

Let us turn to our companion matrix method. Define ideal I = (/i, 72) in the polynomial ring C[z3, Z/^], 
the Grobner basis of ideal I in Lexicographic order Z3 -< Z4 is given by 

GB(/) = (S12S45 + ■S 12 'S 25^^4 — 'Si3S352;4 + S14S45Z4 + 514^252:4, 

S12 + 'S1323 + suZ4 , S45Z3 + S352;4 + 5252324) • ( 4 . 22 ) 


The monomial basis in this Grobner basis is B = { 1 , Z4}. Polynomial reduction of zsB and 2;4B over GB(/) 
gives the companion matrices T^gB = Z3B, T^^B = 2;4B as 


Tz,= 


/ _£i2 

_£M \ 


( 0 

^ ) 

I 513 

513 j 

II 

1 S12545 

514545—513535 1 

[ 512S45 

513535—514545—512525 j 

\ 513S25 

313325 / 


V 314*25 

314325 / 


( 4 . 23 ) 


which are 2x2 matrices, in accordance with the number of solutions of scattering equations. 

We note that the companion matrices actually formally “live” in the ideal I itself by satisfying scat¬ 
tering equations, i.,e.. 


/i —^ • 512 / 2 x 2 + • 513/23 + S 14/24 

—S12 —514 I I [ ^ 

Sl2'S45 ^14S45 —^13'^35 j I _ £l2£45 ■^13^35—■^14^45 —■^12'S25 

and likewise, S43Tz^ + 535^^4 + S2bTz:iTz^ = 02x2- This is, of course, a general property by construction since 
the companion matrices are constructed as multiplication (on a particular basis), so that substituting into 
the defining polynomials would vanish in the quotient ring. The situation is very much analogous to the 
classical result of Gayley-Hamilton that a matrix satisfies its own characteristic polynomial. It is worth to 
emphasize this discussion as 

COROLLARY 2 The companion matrices satisfy the defining polynomials of the given ideal. 

The above corollary shows some kind of equivalence between solutions of scattering equations and com¬ 
panion matrices of monomial basis over the Grobner basis of scattering equations. With these companion 
matrices, we now proceed to compute the trace of the integrands to obtain the amplitude for different 
theories. 


^ = 02x2 , ( 4 . 24 ) 



4.2.1 Scalar Theory 

The 5 -point amplitude of scalar theory is given by 


^5 -Y1 

sol 


z‘^ z^ 
^12^25^51 


1 


1 ^ 


125 

125 


^2 ^2 „2 ^2 „2 
^12^23^34^45^51 


E 


1 

<t|!i(23 - 1 ) = (23 - 


E 

ZS,Z4£SoI 


P{Z3,Z4) 

Q{z3,Z4) 


( 4 . 25 ) 
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where we have used that 

I*|125 _ f '^23 _ ^34 _ / 524 _ ^34 _ 

V (^3 - 1)^ (^3 - Z4y J \ {Z4 - 1)2 {Z3 - 24)2 zj ) {z^ - Z 4 Y 

and as above, defined the appropriate P and Q, which are, explicitly. 


P = Z32(z4-l)^ 

Q = ( 535(^3 - 1)^ + •523^1) {Z3 - Z4f {s4b{z4 “ 1)^ + 524^^1) 


+534 


'S45(^;3 - 1)^^3(^4 - 1)^ 


pz \ (524(^3 ~ 1)^ + '§23(^4 ~ 1)^) + ■535(23 — 1)^(24 — 1)2^ 


(4.26) 


Now, we wish to compute the trace of the matrix PQ~^ upon substituting 23 and 24 by their associated 
companion matrices, instead of snmming over all the complicated solutions of the scattering eqnations. In 
other words, we should replace the variables 23,24 as Tz^,Tz^ in the integrand, i.e., P' = P\z^^Tz^^z 4 ,^t^^i 
Q' = < 3 1 23 ^r 23 ,z 4 ->T ^4 (Hereafter we will always nse P',Q' to denote the matrices after replacing zi to 
Tz^)■ The product of variables 23,24 changes to the product of matrices Tz^,Tz^, and since the companion 
matrices are commutable, their order does not matter in here. Then we should compute the inverse of 
matrix Q', and the final result is given by Tr(P'Q'“^). 

Recalling that the physical poles appearing in the color-ordered amplitude are S12, .§23, .S34, .S45, S15, we 
can define them as the independent Mandelstam variables, and rewrite all the other Mandelstam variables 
in P, Q, T 23 , Tz^ by using following identities: 


■S35 — ■S12 — S34 — ■S45 5 ■524 — Sl 5 — 523 “ ^34 , S25 — 534 “ 5 l5 — S12 , 

514 = 523 — 545 — 515 5 5 i 3 = S45 — S12 — ^523 ■ (4.27) 


After some algebraic manipulation, readily performed by Mathematica, we obtain 


TviP\Tz„Tz,)Q'-\Tz„Tz,)) = 


1 


+ 


1 


+ 


1 


+ 


1 


+ 


1 


(4.28) 


515523 512.534 .515534 512545 523545 

which agrees with the known result [30, 31]. 

Let us further consider an example, corresponding to the two-cycles^ {(1, 2, 3,4, 5), (1, 3, 5, 2,4)}, in 
the language of [30, 31, 36]. Using the CHY-integrand defined by above two-cycles, we have 

1 1 _ Hi(23,24) 


•^5 - ^ 


2 p‘ 

^12^25^51 


sol 


^Il25 ^^12^23^34245^^51 ^132;355;52^24^41 Qi(23,24)’ 


(4.29) 


which is represented by the so-called pentacle diagram (shown in Figure l)from the view of integration 
rules. Using the mapping rule given in [30, 31], the answer is known to be zero. By directly computing the 
trace, we indeed find that Pt{P[Q'{~^) = 0 and confirms this result. 

^Each cycle defines an expression, e.g., Cycle^(l, 3 , 5 , 2 , 4 ) = 1/(213235252224241), and the two-cycles denotes the expression 
given by Cycle^ Cyclej. 
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Figure 1. The pentacle diagram representing the CHY-integrand defined by the two-cycles {(1,2,3,4,5),(1,3,5,2,4)}. 


In fact, for this example, although CHY-integrands of ^5 and are different, after simplification, 
their difference appears only in the numerator, i.e.. 


Pi = - Z 2,){1 - Zi){z2. - Z4) , Qi = Q. 


(4.30) 


Since the trace of matrix is a linear mapping, in particular Tr(Mi -|- M 2 ) = Tr(Mi) -|- Tr(M 2 ), relations 
between results of different integrands should also have hints in the integrand level. For example, let us 
consider the following three CHY-integrands defined by three two-cycles 02 = {(1,2, 3,4, 5), (1,2, 3, 5,4)}, 
03 = 1(1, 2,3,4, 5), (1, 2,4, 5,3)} and 04 = {(1, 2, 3,4, 5), (1, 3, 2, 5,4)}. With some calculations, we find 

P 2 


A(a2) = X] 

sol 

A(a3) = ^ 

sol 

^ 5 ( 04 ) = ^ 


7 ^ ^2 
^12^25^51 

1 

1 

1*1125 

I^ll25 

212223234245251212223235254241 

~2 2 2 
212^25^51 

1 

1 

1*1125 

I^ll25 

2122232342452512122242452:53^31 

P '2 _2 

^12^25^51 

1 

1 


sol 

E 

sol 


E 


sol 


Q 

Q’ 

Q 


sol 1^1125 H22232;345;45^51 2132^32225^54241 

where they share the same denominator Q, but different numerators 

^^2 = 23(2:3 - 2;4)(2;4 - 1)^ , H3 = 23(23 - 1)(2;3 - Z4)(24 - 1 ) , F4 = 2|(23 - 24)(24 - 1 )" 

After putting back the companion matrices, we find that 

1 1 


(4.31) 


(4.32) 


A(«2) = Tr(P'Q'-^) = 
A(«3) = Tr(P'Q'-^) = 


•512^45 

1 


+ 


S 23 S 45 


•512^45 

Realizing that the polynomials have the simple relation 


, A(«4) = Tr(PiQ'-') = 


•523 ■S45 


P^ + Pi — P 2 = Pi , 


(4.33) 
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we obtain the identity amongst these amplitudes as 


TV(P'Q'-i) = Tr((P' + P' - P[)Q'-^) = Tr(P'Q'-i) + TV(P'Q'-i) + Tt{P[Q'-^) 

^ 5 ( 02 ) = ^ 5 (« 3 ) +'^ 5 ( 04 ) + 0 . (4.34) 


Above example demonstrates an idea how to find relations among different amplitudes. Starting from 
different CHY-integrands, we can equalize their denominators by multiplying proper polynomial both 
at the denominator and the numerator. After that, the relations among different amplitudes can be 
understood from the relations among different numerators. 

Let us demonstrate above idea by another example, i.e., the 7-point amplitude-level identity given by 
eq.(3.7) of [31], viz., amplitude obtained from the CHY-integrand 

-i-^- (4.35) 

2^12 ^23 2^34 ^^45 2^56 ^67 2^71 2 2^27 ^74 2^46 ^65 2^53 ^31 

is identical to the sum of following two amplitudes obtained from two CHY-integrand 

-^-^^^-) . (4.36) 

P2^;232:34^;452;56^;672:71 ^;122;56^;372:46 V ^;142:27^;35 ^;252;74^;31 7 

Under gauge-fixing zi = oo, 2:2 = I, - 2^7 = 0 and excluding the 1-st, 2-nd and 7-th scattering equations, the 
Jacobian is 


1 

■ Q ■ ^ ^ 

Thus we can immediately get the numerator of integrand after inserting the above three terms. The first 
term gives 


6 

Hi = 2;32;5(2;4 - 1)(2;5 - 1){zq - 1)(2;3 - zq) Zi{zi - 1) {zi - Zj) , (4.38) 

i=3 3<i<j<6 

while the other two terms give 

6 

P 2 = Z4Z5{Z4 ^ - 1){Z5 - 1){zq - 1){Z3 - Zq) Zi(Zi - 1) (Zi - Zj) , (4.39) 

i=3 3<i<j<6 

6 * 7^5 

Ps = -Z5{Z4 - 1 )( 2:6 - 1)(^3 - ^ 5 )( 2;3 “ Zq) Zi{Zi - 1 ) {Zi - Zj) . ( 4 . 40 ) 

i=3 3<i<j<6 


Note that 


Hi - H 2 - H 3 (4.41) 

6 *7^5 

= 2 ; 5 ( 2;4 - 1 )( 2:6 - 1)(^3 - Z 6 )(Z 4 -Z 5 + Z 3 Z 5 - Z 4 Z 5 ) Zi(Zi - 1) (Zi - Zj) , 

i=3 3<i<j<6 
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while the trace Tr((P{ — P2 — P^)Q' is zero. Note also the following decomposition 


Zi- Z 5 + Z 3 Z 5 - 2:42:5 = (23 - 1)(25 - 24) + 24(23 - 25) , 


SO that we can write Pi — P2 — P3 = P4: + P5, with 


6 

P4, = -Z3{Z3-1){ZI-1){ZQ-1){Z3- ZQ){zi- Z 3 )'^Zi{Zi-l) {zi - Zj) , 

i=3 3<i<j<6 


6 

P 5 = 2425(24 - 1)(26 - 1)(23 - 25)(23 - 25 ) 2i(2i - 1) (zi - Zj) , 

i=3 3<i<j<6 


which correspond to two-cycles 


{( 1 , 2 , 7 , 4 , 6, 5 , 3 ), ( 1 , 2 , 5 , 6, 7 , 3 , 4 )} , {( 1 , 2 , 3 , 4 , 5 ,6, 7 ), ( 1 , 3 , 7 , 2 , 5 , 6, 4 )} 


( 4 . 42 ) 


( 4 . 43 ) 

( 4 . 44 ) 


( 4 . 45 ) 


respectively with Tr(P4Q'“^) = 0, Tr(P5Q'“^) = 0. 

We thus conclude that strictly speaking, the amplitude-level identity between ( 4 . 35 ) and ( 4 . 36 ) is up to 
some CHY-integrands which have vanishing amplitude. More explicitly, the identity ( 4 . 35 )=( 4 . 36 )+( 4 . 45 ) 
holds exactly at the integrand-level, while ( 4 . 45 ) has vanishing final result, so that ( 4 . 35 )=( 4 . 36 ) holds at 
the amplitude-level. This provides the amplitude-level identity an explanation from the basic linearity of 
the trace. 


4.2.2 Yang-Mills theory 

For 5 -point amplitude in Yang-Mills theory, we want to compute 

4YM ^ 'Sp ^ 12 ^ 25^51 Pf ^^ 10 x 10 _ P^'^(23, 24) (AA(\\ 

" h ---23234245251 -^^_^^,Q-(23, 24)’ ^ 

under gauge-fixing 24 = 00, 22 = 1 , 25 = 0 . Let us choose the reduced Pfafhan as 

Pf '^ 10 x 10 = ^^Pf , (4.47) 

21 - 22 

where Til is a 8 x 8 matrix. As in the 4 -point case, the large 24 dependence of Pf ^'I'loxio is I/24, while 
1/(242223234245254) is also l/z\. Together with the factor z‘l2Z23^‘ii i^i numerator, we get a finite integrand 
under the 24 —>■ 00 limit. 

We now follow the standard computation procedure: 

1. Write down the expressions for and Pf'Tioxio, and work out P"^^(23,24), Q^’^(23,24); 

2 . Replace the variables ZiS by companion matrices T^., as P' = P|^._5.'r^., Q' = Q\z^^T:,Z■, 

3 . Compute the inverse of Q' and the trace Tt{P'Q'~^). 
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The result for un-specified helicities is quite lengthy. For illustration, let us consider the 5-point 
amplitude with helicity ^ 5 ^^ 92 ^ dt ^ dt ^ dt)■ polarization vector is dehned as 


^ V2[kr] 


eu(^) = 




(4.48) 


V2[kr] ’ V2{rk) ’ 

and we choose the reference momenta as ri = r 2 = k^,, rz = = k 2 - Thus settled, the only 

surviving products of polarization vectors are e~ {ki)-e^{k 4 ) and e“(A:i) •e'''(A: 5 ). After imposing momentum 
conservation for e^{ki) ■ kj to reduce the ambiguity, we can simplify the 8 x 8 matrix T = as 


0 

fcsfcr 

23-24 

fcafcs 

23-25 

0 

0 

4'16 

£4/03 

24-23 

£ 5^3 

25-23 

24-23 

0 

^ 4^5 

24-25 

eifcr 

21-24 

1 

£. 3^4 

23-24 

3^27 

£ 5 fcl-|-£ 5 fc 3 

25-24 

fcsfcs 

2:5-2:3 

kik ^ 

25-24 

0 

eifc 2 +eifc 4 

21-25 

t 2 ki + e 2 ki 

22-25 

e 3 fcl+e 3 fc 4 

Z3-25 

e 4 fci-l-e 4 fc 3 

24-25 

'I '38 

0 

eifc 4 

21-24 

eik 2 + eik 4 , 

21-2:5 

0 

0 

0 

Cl £4 
21-2:4 

£l £5 

21-25 


€2^4 

22-24 

£3^4 

'^72 

esfcl+es/ca 

25-24 


_ e2ki-\-e2k4 

22- 25 
eski-\-€sk4 

23- 25 
e4/ci+€4fe3 

£4-25 

'I'ss 


(4.49) 


where 


~ ~ eski 63/^4 £ 3^1 + 63/^4 

V 16 = —^61 =-1- 

Zs - Zl Z 3 - Z 4 Z 3 - Z 5 

~ ~ 64/01 , 64/03 64/01 + 64/03 

V 27 = —W 72 =- 1 - 

2:4 - Zl Z4- Z3 Z4 - Z5 

,T, 65/01 65/03 65/01 + 65/03 

V 38 = —^83 =-1- 

2^5 - Zl Z5 - Z3 Z5 - Z 4 


(4.50) 


This greatly simplifies the result of reduced Pfaffian, which reads, after our gauge-fixing, 

T3r/,T, - ^2:4K2,1 +2;46t2,4 - ^t2,l // 

Pf ViOxlO = = -2 -^- -77 -r (2:3 - 2:4)61,5/23,1 K4,1 - (2:3 - 2:4)61,4/23,1 «^5,1 

Dq, Z3Z4(Z4 - 1){Z3 - Z 4 ) 

-(Z3 - ^ 4 ) 61 , 4 / 23 ,1^5,3 - 5;4ei,5K3,4/C4,l + ^^3 61 ,5^3,1^4,3 + ^4e'l,4/«3,4/«5,l) , (4.51) 

where we recall again that = e^ej, Kij = eikj. The factor of scalar part l/( 2 ;i 2 .Z 23 ' 234 - 245 ' 25 i) after gauge 
hxing is , and the Jacobian is the same as in the scalar theory, 


zizi(z3 - 1)^(Z4 - 1)^(Z3 - 2 : 4 )^ 

QYM 


(4.52) 


where is a polynomial of 2 : 3 , 2:4 and Mandelstam variables, and it is also the denominator of integrand. 
The numerator of 1/|4>|J|| cancels the denominator of Pf^T and that of scalar part, leaving a factor 
- 23 ( 2:3 — 1 )( 2;4 — 1) in the numerator. Combined with the numerator Nq, of Pf'TioxiO; they contribute to 
= Z 3 (Z 3 - 1 )(Z 4 - l)Nq,. 
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Then it is straightforward to apply the replacements (T^g, T 24 ) = P™{zz-, (T^g, ) = 

Q'^^{z 3 , Z 4 )\zi^T^., and compute the trace Tr{P '^^To make the computation more efficient, 
we can firstly apply the polynomial reduction of P^^{z 3 , Z 4 ,), Z/i) over GB(/). The remainders 

R{P'^^), are polynomials of Z 4 only, since the monomial basis is { 1 , 2 : 4 }. Then we can proceed by 

replacing 2:4 —)• Tz, for the remainders, and compute the corresponding trace. This gives the same result 
as with the original but the computation would be much faster. With Mathematica, we obtain 


./YM/WYM\-1\ _ r, Cl,5«^2,1^3,1/24,3 + £1,4/22,1/23,4/25,1 — £1,5«^2,1/23,4/24,1 


Tr(P''"^(Q'''’^)“^) = 2 
+ 2 


•512534 

£1,5/22,1/23,1/24,1 + £1,5/22,1/23,1/24,3 — £1,4/22,1/23,1/25,1 — £1,4/22,1/23,1/25,3 

512545 

£1,5/22,4/23,4/24,1 — £1,5/22,4/23,1/24,3 — £1,4/22,4/23,4/25,1 


+ 2 


515534 


(4.53) 


The missing of the pole terms 1 /( 515523 )) 1 /( 523545 ) (terms involving pole S 23 ) is due to the choice 
of polarization vectors. However, the S 23 pole do exist, hiding in / 22 ,i,/ 23 ,i. Directly rewriting the spinor 
brackets for Kij and Sij, and using the Schouten identities we get the famous MHV-amplitude [51, 52] 


Tr(P"“(g'^"")-^) 


(1 2 )" 

(1 2)(2 3)(3 4)(4 5)(5 1) ' 


(4.54) 


4.2.3 Gravity and n-point KLT Relations 

For 5-point amplitude in pure gravity theory, under gauge-fixing zi = 00 , 2:2 = 1) -^5 = 0, we wish to 
compute 


.a? = E 


H2^25^51 


sol 


|4> 


125 

125 


(Pf'4/ioxio)(Pf ^'I'lOxio) = 


^2^25^51 


sol 


|4> 


125 

125 


(Pf'Tioxio)' 


(4.55) 


Let us consider the gravity amplitude ^5 (1 , 2 , 3'’'''', 4++, 5+"''), so that we can use the same reduced 

Pfaffian Pf'Tioxio as in the Yang-Mills case. Here, we do not have the factor of scalar part, but the 
square of the factor of the reduced Pfaffian. The numerator of l/|/h|i 25 cancels the squared denominator 
of reduced Pfaffian z^z^(z 4 — 1)^(.23 — 24 )^, leaving a factor of (zs — 1)^ in the numerator. Hence, we have 
and P° = (zs — with Nq, given in (4.51). 

Thus, all the ingredients have been computed in the Yang-Mills situation above, and we only need to 
work out the trace Tl'(P'^, which gives a lengthy result: 

(1 2 }^ (1 2 f (1 5 ) (3 4 ) [2 1 ]'‘ [3 1 ]^ [4 2 ] [4 3 ]^ [5 1 ]^ -t 971 more termsj 

(1 3 ) (1 4 ) (1 5 ) (2 3 )^ (2 (2 5 )^ (3 4 ) (3 5 ) (4 5 ) [2 1 ] [3 1 ]^ [3 2 f [4 1 ] [4 2 ] [5 1 ] [5 2 ] [5 3 ] [5 4 ] 

where we can see that all poles = 1,..., 5 appearing therein, indicating the colorless structure of 

gravity amplitude. 

This complicated expression can be simplified by non-trivially imposing momentum conservation and 
Schouten identities. Applying the algorithm described in the appendix of [53] , for instance, we can simplify 
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A«(l—,2—,3++,4++,5++) to 


(1 2f [4 3] [5 3] (1 2)® [4 3] [5 4] 

(1 4) (1 5) (2 4) (2 5) (3 4) (3 5) (1 3) (1 5) (2 3) (2 5) (3 4) (4 5) 

, (1 2)^ [5 3] [5 4] 

(1 3) (1 4) (2 3) (2 4) (3 5) (4 5) ' 


(4.56) 


which agrees perfectly with the result given by KLT relation [46-50]. 

More generally, for n-point amplitude, under the usual gauge-fixing zi = oo, Z 2 = 1, Zn = 0, we wish 
to compute 


= 


E 

sol 


^12^2n^nl 


|4> 


12n 

12n 


(Pf'^2nx2n)(Pf'^2n 




x2n) — 


sol 


12 ^ 2 n^nl 
12n 
12n 


| 4 > 


(Pf 


2nx2n) 


(4.57) 


In order to write down the reduced Pfaffian, we need to compute the Pfafhan of a (2n — 2) x (2n — 2) 
matrix, which is quite complicated. Direct computation using the above formula is obviously very difficult, 
just like the direct computation of gravity amplitude by Feynman diagram. So we would like to follow the 
KLT formalism, and compute the gravity amplitude as square of Yang-Mills amplitudes. 

An important property of the reduced Pfaffian is that, it can be expanded [3] as 


Pf'T 


E 

Oi£S„-3 


{zi 


‘S’[«l/3]^n^(l, I3,n,n- 1) _ 

Za2)iZa2 - Zas) ' ' ' iZa„-2 “ ^;n-l)(^;n-l “ Zn)iZn “ Zl) 


(4.58) 


where a, [3 are permutations of labels 2, 3,..., n — 2, and S'fal/I] is the S-kernel. The appearance of 
is a consequence of certain integrand summing over all {n — 2>)\ solutions of scattering equations in the 
original derivation, and in the companion matrix method, it corresponds to the trace of that integrand 
when changing variables to companion matrices. In any event, it is a constant, and can be dragged out of 
the trace. 

Using this expression, we can expand one Pf ^4/ in the gravity amplitude. 


A 


G 

n 


E E 

sol a<^Sn — Z 


/ Piz3,Z4, ■ ■ . , 2 :^- 1 ) \ 

\Q{z3,Z4, . . .,Zn-l)J 


S[a\f3]Al^{l,(3,n,n 

/3eSn-3 


P 

Q 


Z^nZ. 


12 ^ 2 n^nl 


14* 


12n 

12n 


Z1q 2 Za.2Q.3 


Pf ^4^2nx2n 

* —2,77—1^71— 


1 ) , 


(4.59) 

(4.60) 


The trace Tr(P'(T 2 jQ' ^{Tz^)) for the set a gives a, n—1, n), and the summation over permutations 

of a can be taken out of the trace, and we thereby arrive at the KLT relation^. 

'^Note that the ordering of set a(or / 3 ) here defined in [ 3 ] is the reverse of that defined in [ 48 ]. 
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4.3 Six-Point Amplitudes 


We proceed onto six-point amplitudes, i.e., n = 6 . Using the standard gauge-fixing zi = oo, Z 2 = 1, zq = 0, 
Dolan-Goddard’s polynomial form [37] of the scattering equations is given by 

fl = 'S 12 + 513^3 + 514^4 + S 15 Z 5 , (4.61) 

/2 = S123^;3 + 5124^4 + 51252^5 + Sl342:32;4 + Sl352;32:5 + Sl452;42:5 , (4-62) 

/S = S 5 GZ 3 Z 4 546^32^5 + S36242;5 + S26^32;42;5 • (4.63) 

We can thus define the ideal I = (/i,/ 2 ,/ 3 ) in the polynomial ring £.[z 3 ,za,zg\- The degree of ideal I 
is 6 , so according to Bezout’s theorem, it has 6 solutions, though it is not possible to obtain analytic 
expressions for these solutions, as already seen in the 5-point cases. Let us then consider the companion 
matrix method. 

We generate the Grobner basis for I in Lexicographic ordering Z 3 ^ Z 4 ^ z^. Analytically, the explicit 
expression of GB(/) is rather complicated, especially in the presence of so many parameters Sij in the ring. 
By varying the exponents to some high power, the polynomial reduction of the monomials z^^z^^z^^ (with 
Oj from 0 to some finite number, say 20) over GB{I) gives the monomial basis 

B = {l,ZG,zl,zl,zl,zl} . 

The polynomial reduction of 2 ; 3 B, Z 4 B and z^B over GB(/) gives the companion matrices , which 

are 6 X 6 matrices. Again, we need to compute P' = Q' = 

and the final amplitude is given by Aq = Lr{P'Q'~^), without summing over all solutions of scattering 
equations. 

Since the operations we need are multiplication of matrices, taking inverse or trace of matrices, so 
in principle it can be done analytically. However, the symbolic manipulation for n = 6 case is quite 
complicated, especially when taking the inverse of matrix Q' and simplifying the tedious trace result 
in Mathematica, so we introduce random numeric kinematics - i.e., by Monte Carlo assignments of the 
parametres Sij - to get the final result. One will see that, as is customary with coefficient fields in 
polynomial rings, trying a few large prime numbers would suffice very quickly. 


4.3.1 Scalar theory 

We can write the amplitude as 


^6 


^ 12 ^ 26^61 


sol 


I^ll26 


2 2 2 2 2 2 
^ 12 ^ 23 ^ 34 ^ 45 ^ 56^61 


5 - 24)^(24 - 25 ) 32 | 


E 

23,5:4,^5650; 


^( 2 : 34 ^ 4 , 2 : 5 ) 

<3(2:3,2:4,25) 


(4.64) 
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where 


*126 _ 
^126 — 


( $33 $34 $35 

<^43 $44 $45 

V ^53 ^54 ^55 


$ 


126 _ 
126 — 


Det($li) 


(4.65) 


Prime Kinematic Strategy: The idea is the following. Since we know for scalar theory, the final 
result of Aq should take the form 




(4.66) 


where Si- are the independent Mandelstam variables of physical poles si2, S23, S34, S45, sse, sie, S123, S234, 
S345, and the summation is over all possible products of three physical poles, e.g., si2S45^234 ’ 

So in total we have (g) = 84 terms, which we denote as 5i, 1 = 1, 2 ,... 84, and the amplitude is expanded 
as Ae = where Ci is either 0 or 1. 

To each physical pole we now randomly assign a prime number, i.e., we are working with the much 
simpler polynomial ring C\z\ instead of C(s)[ 2 ;]. In this case, the computation of Tv{P'Q'~^) is trivial 
within seconds, and the result as well as 5i’s are all numbers. Next, we shall find the solutions ~ 

for Cl’s. However, doing this by brute-force is impossible since there are 84 Ci’s and each one 
can take 0 or 1, so one would go through all 2®^ possibilities, which is far beyond any computational ability. 

We therefore adopt the following strategy: instead of setting all coefficients to numbers, we can assign 
all physical poles to prime numbers except one pole. For example, we would leave S345, to detect first 
the coefficients of 5i’s which contains the pole Keeping one symbolic variable S345 would extend the 
computation time of Tr(P'Q'“^) up to minutes, but it is still very manageable, while keeping two or more 
symbolic variables would make the computation of Pv{P'Q'~^) in Mathematica very hard for a laptop. 

Let us see the above strategy in action. Setting the kinematics (coefficient variables) as, e.g.. 


S12 = 7 , 523 = 37 , S34 = 79 , S45 = 97 , 

S56 = 131 , S16 = 179 , S123 = 181 , S234 = 223 , 


while leaving S345 free, we get 

Tr(P'Q'-i) 


64909247478 32736 

1878479042622679 ~ 9601739s345 ' 


(4.67) 


Among the Si’s, there are ( 2 ) = 28 terms containing physical pole S345, and the number marked by 
in Tr(P'Q'“^) should be expanded into these 28 terms^. This is thus a problem in Egyptian fractions. By 
going through all 2^® possibilities of Ci, we find the unique expansion 


32736 1 _ / 1 1 1 ^ W 

9601739^ “ V7 x 79 79 x 179 7 x 97 97 x 179/^ ’ 


(4.68) 


5 


In fact, using the compatibility among poles, we can greatly reduce the number of possible combinations of poles. We will 


consider this fact in latter examples. 
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So mapping to the physical poles, we find that 

1 1 


1 


1 


Sl2'S34'S345 Sl6'S34'S345 Si2S45'S345 'S 16 S 45 S 345 

is a part of Aq. 

Now, we try to get more poles. Taking the kinematics as, e.g., 

S 12 = 7 , S 23 = 37 , S 34 = 79 , S 45 = 97 , 

S 56 = 131 , S 16 = 179 , S 123 = 181 , S 345 = 251 , 
while leaving S 234 free, we get 

Tr(P'Q'-^) 


(4.69) 


35960 


13829207594 


68541427S234 302048840293589 

The part marked by the physical pole S 234 can be uniquely expanded as 

35960 1 / 1 1 1 1 \ 1 


+ 


+ 


+ 


thus 


68541427 S 234 V37 x 179 79 x 179 37 x 131 79 x 131/ S 234 

1111 


■Sl6'S23'S234 Sl6'S34S234 S23S56'S234 'S34S56'S234 

is also a part of Aq. With the same procedure, we find that for physical pole S 123 , 

1111 


(4.70) 

(4.71) 

(4.72) 

(4.73) 

Si 2S45'S123 S23'S45'Sl23 Si2S56'Sl23 'S23S56'Sl23 

is also part of Aq. Finally, we need to determine the coefficients Ci of 5i’s without physical poles 
• 5123 , •S 234 , S 345 . There are in total ( 3 ) = 20 terms. Taking the kinematics as, e.g., 

S 12 = 7 , S 23 = 37 , S 34 = 79 , S 45 = 97 , 

S 56 = 131 , S 16 = 179 , S 123 = 181 , S 234 = 223 , S 345 = 251 , 

computing the Tv{P'Q'~^) and extracting the contributions from results (4.69), (4.72), (4.73), the remaining 
result can be uniquely expanded as 

714874 


1 


1 


46539628933 


so the last part for Aq is 


7 X 79 X 131 37 X 97 X 179 

1 1 


•512534556 510523545 

Putting all the above together, we therefore conclude that 

1 1 1 


(4.74) 

(4.75) 


-Tie — — 
+ 

+ 


+ 

512534556 510523545 5i2545Si23 5235455123 5125505123 

11111 

+ 


+ 

+ 


+ 

+ 


1 


+ 

+ 


1 


5235565123 512534S345 5305345345 S125455345 S165455345 

1111 


+ 


+ 


+ 


5165235234 5105345234 5235565234 5345505234 " 

This prime-numeric method can be applied to all the cases of n = 6 amplitudes of scalar theory. 


(4.76) 
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4.3.2 Yang-Mills theory 

For Yang-Mills theory, when n = 6 , we meet the first ’’not so simple” gluon amplitude, i.e., the next-MHV 
amplitude, so it is worthwhile to verify the companion matrix method with this non-trivial example. To 
illustrate, let us consider the split helicity amplitude , §4 , , g^), and choose the reference 

momenta as ri = r 2 = rs = k 4 , r 4 = rs = rg = k^, so that only € 4 ^ 5 , ei^, ? 2 , 5 ) ^ 2,6 are non-zero. The object 
we want to compute is 




YM 

6 


E 


^12^26^61 Pf'^12x12 
l^lll 2^12^;232:34^;452:56^61 


E 

Z 3 ,Z 4 ,Z 5 GsoI 


P'^^{z3,Z4,Z5) 

Q'^^iz3,Z4,Z5) 


(4.77) 


Here both the Jacobian |‘h ||26 reduced Pfaffian Pf ''&i 2 xi 2 are very complicated, so it is almost impos¬ 
sible to compute it analytically. As in the scalar example, we can follow the semi-analytic procedure, 
and set the physical poles as some prime numbers, while keeping 'ki 2 xi 2 (all Qj = eiej, Kij = eikj and kikj 
in 'k) analytic. In this case, the ideal and Grobner basis are just algebraic systems of polynomials with 
integer coefficients, while the elements of companion matrices are rational numbers. So the computation 
is very fast. 

The Jacobian under the chosen gauge-fixing is 

1 ^ zizjzl{z3 - 1)^(Z4 - 1)^(2:5 - 1)^(^:3 " ^4)^(^3 " " Z^^ 

I^ll26 D4.{z3,Z4,Z5) 

where is polynomial in - 2 : 3 , -^ 4 , 2 : 5 . The reduced Pfaffian together with the factor of scalar part give 


^^ (- 2 ^ 3 ? -^ 4 ? ^5 ? ^i,j 5 ^ij ? kikj ) 

2:32:42:5 (^^3 - 1 )^( 2;4 - 1)(25 - 1)(2;3 “ 24 )^(^3 “ 25) (2:4 - 2:5)^ 
under the chosen gauge-fixing for some polynomial numerator So we have 


(4.79) 


P™( 2 : 3 , 2:4, 25 ) = 2:32:4(24 - 1)(25 - 1 )( 2:3 - 2 : 5 )^^ , (5™(23, Z4, 25 ) = 7?$ . (4.80) 


Note that N 4 , originates from the Pfaffian of a 10 x 10 antisymmetric matrix, where by definition, each term 
in the Pfaffian is a product of five elements in the matrix. So each term in Nq, is a product of five elements 
selected from?jj, kij, kikj, combined with a monomial of 2 : 3 , 2 : 452 : 5 , for example, 2232 |? 2 , 6 ^ 1 , 2 ^ 3 , 4 / 24 , 1 ^ 5 , 1 . 
Finally we can take the replacement p''^^ = and compute the trace 

Tr(P'YM(Q/YM)-l)_ gjyg ]3 g^g 

1 44 1 

Tr(T’'YM(g'YM^ ) = 0 X 41149 ^2,6«:i,2/23,2/24, 1 ^ 5,1 - 5/21,3/23,2/24,1/26,2 

-|- 500 more terms . (4.81) 


Using the techniques shown in scalar (pt' theory, we can uniquely decompose the rational numbers as 


44 


1 


+ 


1 


1 


1 


6141149 7 X 181 X 131 181 x 37 x 131 


877307 181 X 37 X 131 


(4.82) 
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so we can conclude that 


Tr(P'''M(Q'™) — ( -1-) 62 ,6^1,2^3,2«^4,1^5,1 

^Sl 2 Sl 23 'S 56 'Sl 23 'S 23 'S 56 ^ 

-e2,56vi,3K3,2«^4,i6v6,2 + 500 more terms . ( 4 . 83 ) 

'Sl 23 'S 23 'S 56 

Rewriting them as spinor products and applying the simplification algorithm for spinor expression, we get 
a one-page long result, which remarkably agrees with the known answers [ 54 , 55 ]. 


4.4 Seven-Point Amplitudes 


The companion matrices T^. are simultaneously diagonalizable, and according to Stickelberger’s theorem, 
the complex roots Zi of ideal I are the vectors of simultaneous eigenvalues of the companion matrices T^.. 
Thus when they are evaluated in the diagonal form, the matrices P' = , Q' = , P'Q'~^ 

are also diagonal, and it builds the one-to-one mapping between diagonal elements of (n — 3 )! x (n — 3 )! 
matrix P'Q'~^ and the integrand P/Q evaluated at the (n — 3 )! complex solutions of scattering equations. 
To demonstrate this, let us go through a 7 -point example of scalar cf)^ theory. 

As usual, let us gauge fixing = oo, 22 = 1 , 2:7 = 0 , and the amplitude is given by 


^7 = Y^ 

sol 


~2 2 2 
^ 12 ^ 27^71 


1 


14* 


127 

127 


z‘^ z‘^ z‘^ z‘^ z‘^ z‘^ z‘^ 
^ 12 ^ 23 ^ 34 ^ 45 ^ 56 ^ 67^71 


E 

2 : 3 , 2 : 4 , 25 , 


P{z 3 ,Z 4 ,Z 5 ,Zg) 

Q(z 3 ,Z 4 ,Z 5 ,Ze) 


( 4 . 84 ) 


where 


<j,127 

^127 


/ 4>33 <I >34 4>35 4>36 \ 

4>43 <I>44 4>45 4>46 

4>53 4>54 4>55 4>56 

\ 4)03 4>64 ‘hes 4>66 y 




127', 

I27J 


( 4 . 85 ) 


The Dolan-Goddard polynomial form [ 37 ] of the scattering equations is given by 


fl = S12 + S13Z3 -h 5442:4 -F 5452:5 -h 546^6 , (4.86) 

/2 = 54232:3 + 51242:4 + 542525 -I- 542020 

+ 54342324 -F 5435^32:5 -h 54302:32:6 -b 54452:4^5 -b 54402:42:0 -b 5450^52:0 , (4.87) 

/3 = 542342324 -b 542352325 -b 542362320 “b 542452425 “b 542462420 “b 542562020 

+ 54345232425 -b 54346232420 -b 5435623202:0 -b 5445624202:0 , (4.88) 

/4 = 5572:32:42:5 -b 5572:32:42:6 -b 5472:32:52:5 -b 5372:4^52:0 -b 5272324252:0 . (4.89) 


We can define the ideal I = {fi, f2, f3, f4) in polynomial ring C[2:3, 2:4, 2:5, 2:0], and generate the Grobner 
basis of I in Lexicographic order 2:3 -< 2:4 -< 2:5 -< 2:0. The degree of ideal I is 24 , so the variety of ideal I is 
given by 24 point solutions for which there are no closed form solutions. 
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Let us set the kinematics (all physical poles) as some prime numbers, 

•512 = 5 , S23 = 37 , S 34 = 43 , S 45 = 61 , S56 = 97 , se? = 101 , S 17 = 139 , 

•5123 = 151 , S 234 = 163 , S 345 = 191 , S456 = 211 , §567 = 223 , 5671 = 251 , 5712 = 263 (4.90) 

in the following computation. The solutions of scattering equations fi = 0,i = 1, 2, 3,4 requires computing 
the roots of equations of degree 24, which has no closed form in radicals. Doing it numerically, we get 24 
solutions 

soil : 23 = 20.9071 , 24 = 1.66835 , 25 = 7.08198 , 25 = -64.2332, 

S0I2 : 23 = 1.4223 - 0.318993i , 24 = 12.204 - 5.48743i , 25 = 0.342956 - 0.477119i , 23 = 51.9097 - 32.886i, 

soh ■■ Z3 = 1.4223 + 0.318993i , 24 = 12.204 + 5.48743i , 25 = 0.342956 + 0.477119i , 23 = 51.9097 + 32.886i, 

5014 : 23 = 27.2316 , 24 = 1.76178 , 25 = 13.0497 , 23 = -12.5157, 

5015 : 23 = 1.34598 , 24 = -3.76733 , 25 = -1.28282 , 23 = -56.7763, 

soZs : 23 = 4.92534+ 1.82303i , 24 = 2.04236 + 0.470521, 25 = 0.12331 + 0.733661, 23 =-36.88 - 1.748571, 

50^7:23 = 4.92534- 1.823031, 24 = 2.04236- 0.470521, 25 = 0.12331 - 0.733661, 25 =-36.88 + 1.748571, 

sols : 23 = -11.2804 , 24 = 3.5116 , 25 = -6.80042 , 25 = -1.26394, 
soIq : 23 = 1.19261 , 24 = -8.20104 , 25 = 3.07784 , ze = 6.22689, 

soho : 23 = 1.18325 + 1.937451 , 24 = 0.29585 - 0.486391 , 25 = 0.56997 + 1.110081 , 25 = 0.15405 - 0.393591, 

so/11 : 23 = 1.18325 - 1.937451 , 24 = 0.29585 + 0.486391 , 25 = 0.56997 - 1.110081 , 25 = 0.15405 + 0.393591, 

so/12 : 23 = -4.76521 , 24 = -3.05026 , 25 = -1.6908 , 25 = -0.488528, 
so/13 : 23 = 0.576445 , 24 = -3.05135 , 25 = 1.14498 , 25 = 0.712806, 

so/14 : 23 = 1.78095 + 0.416391 , 24 = 2.1103- 0.606631, 25 = 0.52752 + 0.299271, 26 = 2.3283 - 1.390611, 

so/15 : 23 = 1.78095 - 0.416391 , 24 = 2.1103 + 0.606631, 25 = 0.52752 - 0.299271, 26 = 2.3283 + 1.390611, 

S0/16 : 23 = 1.86192 , 24 = 0.877999 , 25 = 0.795994 , 26 = 0.601979, 

so/17 : 23 = 1.65547 , 24 = 1.9848 , 25 = 0.493798 , 26 = 2.31186, 

S0/18 : 23 = 0.327576 , 24 = -0.0855936 , 25 = 0.212916 , 26 = -0.0559545, 
so/19 : 23 = 0.307828 , 24 = 0.645287 , 25 = 0.0420044 , 23 = 0.46483, 

so/20 : 23 = 0.174313 , 24 = 0.120642 , 25 = 0.0855984 , 26 = 0.0445606, 

so/21 : 23 = 0.031819 , 24 = 0.15022 , 25 = 0.00455376 , 23 = 0.0545382, 
so/22 : 23 = 0.0191033 , 24 = 0.0765079 , 25 = 0.0145344 , 26 = 0.00921803, 

so/23 : 23 = -0.100486 , 24 = -0.0950558 , 25 = -0.00857861 , 26 = -0.0892275, 

so/24 : 23 = -0.0162083 , 24 = 0.0167369 , 25 = 0.00970032 , 26 = -0.0265855 , 
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Let us now turn to the companion matrix method. The monomial basis over GB(/) is given by 24 
elements 


T 3 ^ ri _ .2 3 4 5 6 7 8 9 10 11 12 


US 


.,14 


U5 .,16 


,17 




,18 „19 

Zr 


,20 


,21 


,22 


,23 


} 


(4.92) 


Accordingly, by polynomial reduction of ZiJi,i = 3,4, 5,6 over GB(/), we can get the companion matrices 
Tzi,i = 3,4, 5,6, which are 24 x 24 matrix and satisfying Tz^Ji = z^B. In order to compute Aj, we can 
proceed as usual by computing TT{P\z^^Tz^Q\2^-^Tz )’ result is 


Tr(P'g'-i) 


19260317055974762778118 

9649229470008137021319652355 


1.99605 X 10“® , 


(4.93) 


which agrees with the numeric result given by summing over all solutions of scattering equations. Again 
we see that, since the computation only involves basic manipulations on matrix, we are able to get the 
closed form result, and show that the final result is rational functions of Mandelstam variables. 

The companion matrices are simultaneously diagonalizable. We can choose T^g and compute its eigen¬ 
vectors, since T^g is the simplest companion matrix by definition. Such computation involves finding the 
roots of equations of degree 24, which prohibits analytic solution. Now, T^g has 24 column eigenvectors 
UiA = 1) • • •) 24, and from them we can define the transformation matrix U = {ui, ..., W 24 ) 24 x 24 - Then 
T^. = U~^Tz^U, 7 = 3,4,5, 6 are all diagonal matrices, explicitly given as 

diagfT^) = {20.9071 , 1.42230-h 0.31899i , 1.42230 - 0.318991, 1.345982, 4.92534- 1.823031, 

4.92534-^ 1.823031, 27.2316, 1.192609, 1.78095 - 0.416391, 1.78095-^0.416391, 

1.65547, - 11.28035, 0.576445, 1.86192, -4.76521, 0.307828, 

1.18325 - 1.937451 , 1.18325 -h 1.937451 , - 0.1004864 , 0.327576 , 0.0318190 , 

0.174313 , - 0.0162082 , 0.0191033} , 


diag(T{() = {1.66835 , 12.20402 -^ 5.487431, 12.20402 - 5.487431, - 3.76733, 2.04236 - 0.470521, 

2.04236-^0.470521, 1.76178, - 8.20104, 2.11030-^0.606631, 2.11030- 0.606631, 

1.98480 , 3.51160 , - 3.05135 , 0.877999 , - 3.05026 , 0.645287 , 

0.295854-^0.4863861, 0.295854- 0.4863861, - 0.0950558, - 0.0855936, 0.150220, 

0.1206420 , 0.0167369 , 0.0765079} , 

diag(T{() = {7.08198 , 0.3429560.4771191 , 0.342956- 0.4771191 , - 1.282815 , 0.123310 - 0.7336581 , 

0.123310-^0.7336581, 13.04970, 3.07784, 0.527517- 0.2992731, 0.527517-^0.2992731, 
0.493798 , - 6.80042 , 1.144984 , 0.795994 , - 1.69080 , 0.0420044 , 0.569967- 1.1100811 , 

0.569967-^ 1.1100811, - 0.00857861, 0.212916, 0.00455376, 0.0855984, 0.00970033, 

0.0145344} , 
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and 

diag(T4) = {-64.2332 , 51.9097 + 32.8860i , 51.9097 - 32.8860i , - 56.7763 , - 36.8800 + 1.7486i , 

-36.8800 - 1.74861, - 12.51570, 6.22689, 2.32830 + 1.390611, 2.32830- 1.390611, 

2.31186, - 1.263937, 0.712806, 0.601979, - 0.488528, 0.464830, 0.154052 + 0.3935881, 

0.154052 - 0.3935881, - 0.0892275, - 0.0559545, 0.0545382, 0.0445606, - 0.0265855, 

0.00921802} . 


It can be checked directly that, each set of diagonal elements {Tz(^)i^i} corresponds 

to a set of solution of scattering equations. Thus each diagonal element of matrix 

P'{T^.)Q'~^{T^.) is identical to the integrand P/Q evaluated at one solution of scattering equations, and 
the equivalence between results of these two methods is obvious. 

With the arithmetic result, it is possible to determine the terms appearing in amplitude by setting 
appropriate kinematics. In fact, in this example, we know that the result should be the sum 

^7 = E- 

... 'SllSl2'Sl3Sl4 


where naively the summation is over all possible products of 4 physical poles s,., i.e., = 1001 terms 

with Cl, 1 = 1,... 1001, being either zero or one. By choosing 1001 different group of kinematics for physical 
poles, we get 1001 linear equations of (4.94), and solving them gives the Cj. 

Indeed, the number of terms grows very fast with n in (4.94). The number of independent poles is 
_ (n-iXn- 2 ) _ ^ massless theory, while the number of possible terms in the expansion is ■ For 
n = 8, the number is 15504, and for n = 9 the number is 296010. So it is not very doable when n is 
large. However, for theory, the number of color-ordered diagram is much smaller, and the counting is 
—So for n = 7, the possible terms appearing in (4.94) is 42 (an auspicious number). For n = 8, 
the number is 132, and for n = 9, the number is 429, etc. If we restrict to the 42 possible terms in (4.94), 
then it is possible to compute the coefficients Ci by choosing one set of kinematics. 

One can let each physical pole be assigned a random prime number, and compute Tr{P'Q'~^) and then 


let Mathematica go through all 2^^ possibilities of Ci’s to find the summation 


= Tr(P'Q 




If the prime numbers in kinematic variables are distributed randomly in a very large scale, e.g., primes 
between 2 to 10000, then usually we can find one unique solution for c, in the spirit of Egyptian fractions. 
This enables us to do one computation and fix all coefficients. 

For example, let us compute 


•^7 - X] 


■P P- -P 
^ 12 ^ 27^71 


1 


sol 


|<T)|127 

I^ll27 


2^12^^23^34245 2^56^^67^712122:245^45^572762:63^^31 


(4.95) 


With the kinematics shown in (4.90), we find the unique decomposition 


Tr(P'Q'-^) = 


284 


1 


+ 


1 


1037296765 5 x 61 x 101 x 151 5 x 101 x 151 x 233 ’ 


(4.96) 
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which indicates that 


A'j — I , (4.97) 

Sl2'S45'S67Sl23 •Sl2'S67Sl23'S567 

and agrees with the result given by CHY mapping rules [30, 31]. 

There is a way to directly determine whether a certain term - ^ - is present in the result or not, 

S,^Si2Si3Si4 

by setting the kinematics = a, Si 2 = a^, 5,3 = a^, 5,4 = a®, and others random primes not equaling to 
a. If this term exists, then the denominator has a factor . Again in the A!rj example, if we instead set 
S 12 = 5, S 45 = 5^, S67 = 5^, S 123 = 5®, then the result is ? thus ^ ^ \ ^ is a term in A'j. However, if 

o Xzzo 1 S 12 ^13 ^14 

we set S 12 = 5, 555 = 5^, S 67 = 5^^, S 123 = 5®, then the result is 5i3xmt<223 - indicates that 34233^^676123 
is not a term in A!j, while the 5^^ factor indicates that possible terms involving 342^676123 exist, which 
provides further information for detecting other existing terms. By this way, we can check all possible 
terms by setting kinematics for each one. 

The number of solutions for scattering equations grows as (n — 3)!, while the companion matrix grows 
as (n — 3)! x (n — 3)!. When n = 8, we need to invert the matrix Q 120 X 120 ’ at n = 9, the matrix 

^720x720’ This sets the limitation on the computation of higher n. 

5 Conclusions and Outlook 

In this paper, motivated by the explanation of equivalence of different integrands in the CHY setup, we 
propose a new method using companion matrices, borrowed from the study of zero-dimensional ideals 
in computational algebraic geometry, to evaluate the integrand. One advantage of the method is that 
the rationality of final integral is obvious. Thus although our method may not be as efficient as the one 
proposed in [30, 31, 36], it does give a new angle to study the important problem of scattering amplitudes. 

As shown in the plethora of examples, when the number of external legs grows, the analytic expression 
of companion matrix becomes harder. In fact, when n > 6 , the best way to do it is by assigning the 
kinematic variables to random prime numbers in order to reconstruct the analytic result. The salient 
feature of our method is that it is purely linear-algebraic, involving nothing more than finding the inverse 
and trace of matrices. The linearity of the trace, for example, was demonstrated to immediately lead to 
non-trivial identities in the amplitudes. 

Now, since the physical problem is very symmetric as can be seen by the polynomials given in (2.15) 
and (2.16), one is confronted with an immediate mathematical challenge. If we could analytically find, say 
by induction, the Grobner basis and subsequent monomial basis for the polynomial form of the scattering 
equations in some appropriate lexicographic ordering, then one would find a recursive way to construct the 
companion matrix explicitly, much like the recursive construction of tree-level amplitude by using BCFW 
deformation [16, 17]. Working out this construction is hard but worthwhile, as it would give explicit 
analytic results for the amplitudes and provide a deeper understanding of the CHY formalism. 
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